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Effect of disorder studied with ferromagnetic resonance for arrays of 
tangentially magnetized sub-micron Permalloy discs fabricated by nanosphere 
lithography 

N. Ross/'B M. Kostylev/'E3and R. L. Stamps^ 

School of Physics, University of Western Australia, Crawley, WA, Australia 

(Dated: 5 July 2010) 

Tangentially magnetized trigonal arrays of sub-micron Permalloy discs are characterized with ferromagnetic 
resonance to determine the possible contributions to frequency and linewidth from array disorder. Each array 
is fabricated by a water-surface self-assembly lithographic technique, and consists of a large trigonal array of 
700 nm diameter magnetic discs. Each array is characterized by a different degree of ordering. Two modes are 
present in the ferromagnetic resonance spectra: a large amplitude, 'fundamental' mode and a lower amplitude 
mode at higher field. Angular dependence of the resonance field in a very well ordered array is found to be 
negligible for both modes. The relationship between resonance frequency and applied magnetic field is found 
to be uncorrelated with array disorder. Linewidth is found to increase with increasing array disorder. 

Keywords: ferromagnetic-resonance, linewidth, nanosphere-lithography, sub-micron-disc 
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I. INTRODUCTION 

There is currently an intense interest in the magnetic 
properties of nanoscale and sub-micron discs of low as- 
pect ratio. Such discs have potential applications in data 
storageji"— spintronics;^ and medicinej^i^ and are other- 
wise viewed as simple model systems by which the prop- 
erties of more exotic nanostructures can be understood. 
In many fundamental studies and in data storage and 
spintronics applications in particular, the dynamic mag- 
netic properties of the disc structures are particularly im- 
portant. Where arrays of discs are densely packed, dipole 
coupling between discs can have a significant impact on 
dynamic behaviour. There have been a variety of exper- 
imental and theoretical studies of the dynamic magnetic 
properties of arrays of dipole-coupled magnetic discs, 
with recent publications concerning resonance frequency 
position with respect to applied magnetic field, ^''"^'^ mode 
structure in the vortex state^^ii^ and mode structure 
for in-planei^ii^"— and out-of-planeii'— ^ applied mag- 
netic fields. Recognition of the importance of spin- wave 
damping for many applications has led to a number of 
studies into ferromagnetic resonance (FMR) linewidth 
broadening in such systemSj^^— and recently in a study 
of microwave- assisted switching of magnetic array nano- 
elementSjSi To date, however, an experimental study of 
the effects of array packing on damping of spin wave 
modes observable with FMR has been lacking. 

One approach to such a study is to measure FMR 
linewidth damping in arrays that were patterned by a 
'top-down' method like Focussed Ion Beam lithography 
in such a way that they had identical disc geometries but 
varied pitch (inter-disc spacing). Such fabrication tech- 
niques are time-intensive and therefore often only suit- 
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able for the production of small arrays (< 1 mm^). The 
FMR responses of small arrays are correspondingly small. 
Such small signals are often beyond the sensitivity of con- 
ventional FMR techniques, so that the previous experi- 
mental studies of damping in arrays of nano-discs uti- 
lized time-resolved Kerr microscopy^Si^ and 'meander- 
line' FMR^^ for the detection of spin wave modes . In 
this study, nanosphere lithography has been used to fa- 
cilitate the production of large-area arrays, thus allowing 
measuring using conventional Vector Network Analyzer 
FMR (VNA-FMR). 

'Nanosphere lithography' refers broadly to techniques 
in which the tendency of colloidal particles of sub-micron 
or nanoscale size to self-assemble on hydrophilic surfaces 
is exploited to produce a lithographic mask; the rele- 
vant literature is replete with various examples of such 
techniques fSi^iS The most pronounced limitation of such 
techniques is the difficulty of achieving long-range order- 
ing in the colloidal mask. 

In this study, this tendency towards partial disorder 
was used to fabricate a series of array samples with var- 
ied degrees of array ordering. A nanosphere lithographic 
technique involving assembly of colloidal particles into a 
monolayer on a water surface^ was used to produce ar- 
rays of sub-micron magnetic discs with deliberately var- 
ied degrees of array order but the same local disc-to-disc 
spacing and disc geometry. The FMR mode structures 
and linewidths of this series of array samples were mea- 
sured by VNA-FMR. 



II. EXPERIMENT 

Four sub-micron magnetic disc array samples with dif- 
ferent degrees of ordering were fabricated by a technique 
derived from that of Weekes et al.'^^ Three continuous 
NigiFeig (Permalloy) films, thickness L — 27 ± 3 nm, 
were radio-frequency magnetron sputtered on silicon sub- 
strates with 40 nm Ta seed layers: these films are de- 



noted fl, f2, and f3 in this study. Sputtered thickness 
was verified by the measurement of step edges of inde- 
pendently sputtered films using white-light optical pro- 
filometry. Each of these films was placed on a raised 
platform in a narrow bath filled with 18 Aff2 • cm wa- 
ter, so that the Permalloy face of the film was under 
the water surface but as close to it as possible. Sodium 
Dodecyl Sulphate (SDS) was added to this bath to a con- 
centration of 150 fig L~^: it has been demonstrated in 
the literature that SDS can aid the assembly of reagent 
microspheres on a water surface4^i^ 780 nm diame- 
ter carboxylate-terminated polystyrene reagent micro- 
spheres obtained from Duke Scientific were mixed with 
reagent-grade ethanol into a 1:2 microsphere solution : 
ethanol mixture by volume. This mixture was intro- 
duced to the water surface via a hydrophilic glass slide 
previously cleaned in SDS. Tefion sliders placed laterally 
across the water bath were used to gently agitate the 
water surface, aiding in the self-assembly of well ordered 
trigonal monolayers of microspheres. These sliders were 
used to position the monolayers above the film, and the 
water drained quickly from the bath in order to transfer 
the monolayer to the film surface. Varying the degrees 
of agitation and of monolayer compression during the 
transfer stage allowed control over the degree of order- 
ing of the resulting lithographic mask. The assembly of 
the lithographic mask on the water surface is depicted 
diagrammatically in Figure [TJ 
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FIG. 1. Schematic cross-section representation of the depo- 
sition of the lithographic mask, followed by agitation of the 
water surface with laterally placed teflon sliders to aid assem- 
bly of trigonal monolayers. 



Once masked with the trigonal array of microspheres, 
several 5x8 mm^ sections were cut from each of the par- 
ent films fl, f2, and f3: one such section from each parent 
film was left masked but was not patterned; the other sec- 
tions were patterned. In this study, un-patterned sections 



of the parent film are denoted flc, f2c, and f3c. Four films 
with different degrees of mask ordering were patterned: 
in this study they are denoted fla, f2a, fSa, and f3b to 
make clear the identity of the film from which each was 
patterned. Each film was patterned by placing it in the 
target position of a radio-frequency sputtering chamber 
which was evacuated to a base pressure of 1.0 x 10~^ Torr. 
Samples were reactively ion etched in an oxygen atmo- 
sphere of 100 mTorr at a power density of 0.5 W cm~^ in 
order to reduce the diameter of the mask units. Follow- 
ing etching, each sample was argon milled at a pressure 
of 25 mTorr and power density of 4.25 W cm~^ to remove 
Permalloy not masked by microsphere material. 




FIG. 2. AFM image representative of the microscale character 
of the disc arrays; the high contrast regions in the centre of 
the discs indicate that some of the polystyrene caps remain. 
This particular image is a 5 x 5 /xm image of sample f3b. 



Characterization of disc geometry was achieved with 
tapping-mode Atomic Force Microscopy (AFM): an ex- 
ample AFM image is shown in Figure [2j With the scan 
parameters used, the AFM could be used to resolve struc- 
tures down to 20 nm. The high structures in the cen- 
tre of the discs in the AFM image are the remains of 
polystyrene caps. The discs in all of the arrays have very 
smooth edges and the local order is very good. Values for 
average diameters were extracted by taking cross-sections 
of rows of dots from several such images from each sam- 
ple. 

Static magnetic properties of the array samples were 
measured using a SQuID magnetometer. Hysteresis 
loops were recorded at T = 295 K. Each sample dis- 
played the double 'closed loop' characteristic of vortex 
nucleation, movement, and annihilation in sub-micron 
magnetic discs. Shown in Figure [3] is the loop corre- 



0.0004 



_ 0.0002 

i 0.0000 



-0.0002 



-0.0004 







-2000 -1000 



2000 




FIG. 3. SQuID-measured hysteresis loop of f3b, typical of 
the loops recorded for all the patterned films; the fit lines in 
this figure are linear interpolations that serve as a guide for 
the eye. The hysteresis loop shows the two closed loops char- 
acteristic of vortex nucleation, movement, and annihilation. 
The narrowness of the 'pinching off' of the two loops provides 
confirmation for the whole sample of smooth dot boundaries 
seen in Figure (2] 
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spending to sample f3b, the same sample shown in the 
AFM image in Figure [2J Hysteresis loops of the other 
samples displayed the same characteristics. 

Array regularity of each patterned sample was mea- 
sured using Scanning Electron Microscopy (SEM). Four 
evenly spaced 3600 x magnification images were taken 
from each 1 mm section of the centre line of the long 
axis of the film. Each of these sets of four images were 
stitched together to make a composite image, and the 2- 
dimensional Fourier transform of each composite calcu- 
lated. An example of such a single SEM image with the 
Fourier transform of the composite image from a differ- 
ent section inset is shown in Figure |H The length of each 
sample was L = 8 mm. The Fourier transform from each 
1 mm section of the sample was used to measure A(j)i , the 
variation of lattice angle over that one millimeter. The 
eight separate A0i values were averaged to yield (f>\ the 
average variation in lattice angle per unit millimeter of 
the film: (j)' = X]i=i ^4>ilL. </>' was used to characterize 
the array ordering of each patterned sample. 

Each of the continuous and patterned films was charac- 
terized using microstrip- waveguide VNA-FMR. The sam- 
ple was placed face-down on an 8 milli-inch copper mi- 
crostrip waveguide in an in-plane saturating magnetic 
field. A network analyzer provided microwave excita- 
tion to the waveguide and measurement of transmission 
parameter 5*21. Negligible reflections allowed 5ii to be 
ignored."''* In contrast to how VNA-FMR is typically per- 
formed, the excitation frequency / was fixed and the ap- 
plied magnetic field H varied, in analogy to a cavity FMR 
measurement. The measurements were performed at 1 
GHz intervals in the domain 7-20 GHz. A typical spec- 
trum is shown in Figure [S) Two modes can be resolved: 
a large amplitude mode, and a smaller amplitude mode 
at higher field (lower frequency). Resonance fields H^cs 



FIG. 4. An example SEM image with inset 2-D Fourier trans- 
form of a montage of 4 evenly spaced images taken over a one 
millimeter section of the sample: these particular images are 
also of sample fSb. Note that the SEM image shown was not 
among the images used to produce the composite EFT. Im- 
ages similar to the inset from each 1 mm section of the sample 
were used to calculate the average variation in lattice angle 
per unit length, 0', for each sample. 



and linewidths AiJ were extracted from these measure- 
ments by least-squares fitting Lorentzian curves to the 
troughs in S'21, using at least twenty-five data points on 
each side of the minimum. / vs i?res data for continuous 
samples fl-3c was fitted with the Kittel equation^ with 
demagnetizing factors N^ = Nz = 0, Ny = Att: 



u^l-i" ={H + Hk + {Ny - N,)Ms)x 
{H + Hk + {N^-N,)Ms). 



(1) 



and saturation magnetizations AttMs extracted. Cav- 
ity FMR measurements were performed to verify the ob- 
served mode structure. 

After characterization of all samples, sample f3b was 
cut to obtain a smaller sample (5x5 mm^) with a much 
higher degree of ordering. This sample, denoted f3b*, 
was then characterized in the same way as the other 
four samples. Additionally, FMR spectra were recorded 
for this sample at an excitation frequency of 10 GHz 
for various relative angles (j) between the static applied 
magnetic field and the lattice vector defining the line 
between nearest neighbours. In addition to providing 
a more ordered sample, comparison of FMR results be- 
tween f3b and f3b* allowed the possibility of parent film- 
or patterning-caused differences between array samples 
to be more reliably discounted. 
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FIG. 5. Plot of normalized intensity of transmission ampli- 
tude 521 vs applied magnetic field H for sample f2a at an 
excitation frequency of 16 GHz. In addition to the large- 
amplitude fundamental mode, there is a small mode at higher 
applied field. 



III. RESULTS AND DISCUSSION 



A. Film structure 



The disc diameters as measured by AFM, array order- 
ing values as measured by SEM, and parent film magne- 
tization values as measured by VNA-FMR are presented 
in Table HI The variation in disc diameter is on the or- 
der of five percent, which is the manufacturer's quoted 
uncertainty on the size of the reagent microspheres used 
and therefore the best uniformity achievable by this fab- 
rication technique: the disc sizes of all four samples were 
identical within the limits of experimental uncertainty. 
The narrowness of the 'waist' of the SQuID hysteresis 
loops of the samples (Figure ^ showed that there were 
very few uncompensated magnetic moments at rema- 
nence, likely indicating very little roughness to the disc 
boundaries. 

Each of the five patterned samples was distinguished 
by a unique ordering parameter </)', although there was 
some overlap between f3a and f2a at the limits of exper- 
imental uncertainty. Each of the three parent films' sat- 
uration magnetizations Ms were slightly different. The 
linewidths of these films at a given frequency, /S.H were 
also different. For this reason, VNA-FMR data for the 
parent continuous films f3, f3, and fl are included in the 
tables in this section. In Figure [5] the linewidths for the 
parent continuous films are presented alongside those of 
the corresponding patterned films. In the comparison of 
inhomogeneous broadenings between films shown in Fig- 
ure [5] these parent film linewidths have been accounted 
for. 



Film 


d [nm] 


<^'[° mm-^ 


AtxMs [kOe] 


f3c 


- 


- 


8.49 


f3b* 


695 ± 28 


2.6 ±0.4 


- 


f3b 


695 ± 28 


6.0 ±0.8 


- 


f3a 


703 ± 37 


9.4 ± 1.1 


- 


f2c 


- 


- 


8.69 


f2a 


697 ± 31 


11.3 ±1.7 


- 


flc 


- 


- 


8.85 


fla 


699 ± 28 


19.9 ±2.1 


- 



TABLE I. Table showing average disc diameter, d, array vari- 
ation per unit length </>', and saturation magnetization AnMs 
for the samples used in this study. 



B. FMR mode structure 

A similar mode to the high-field mode in Figure 
[5] has been calculated to exist for isolated discs of 
larger aspect ratio L/d that those fabricated for this 
study, using OOMMF'^" and linearized micromagnetic 
approaches j22i^ In these publications the mode is calcu- 
lated to be located in the 'end' of the disc as defined by 
the direction external magnetic field. These two studies 
also report experimental measurements of these modes 
in closely packed disc systems. Such modes have also 
been measured elsewhere in systems of comparable disc 
dimensions to those fabicated for this studyJ^ii^ VNA- 
FMR lacks the kind of spatial resolution of techniques like 
magnetic resonance force microscopy, and in this study 
the possible confinement of the high-field mode to specific 
regions of the discs could not be investigated directly. 

Several possibilities for the origin of the mode were 
ruled out experimentally. The possibility that the mode 
was caused by a non-uniform excitation field result- 
ing from microwave screening by eddy currents in the 
sample^^ was ruled out by measurement of the mode 
structure by cavity FMR. In an FMR cavity, the excita- 
tion field is very uniform. Figure IHl shows the cavity FMR 
spectrum. The higher field mode is still present, indicat- 
ing that it does not result from non-uniform microwave 
excitation. The amplitude of this mode as measured by 
VNA-FMR did not vary significantly between samples of 
different degrees of ordering. This ruled out the possibil- 
ity of the high-field mode being the result of a collective 
excitation occurring at discontinuities or around vacan- 
cies in the array, which are more numerous in more disor- 
dered samples. Finally, the out-of-plane saturated FMR 
spectrum-not included here for brevity-showed up to five 
well resolved modes in the structure expected for cylin- 
drical symmetry^ suggesting that the high-field mode 
was not the product of some repeated non-cylindrical ele- 
ment. Remaining possible explanations for the high- field 
mode are that it is an 'end' or 'edge' mode as in Refer- 
ences |T3,|2J, andlSO. or that it arises out of or is modified 
by some collective effect of the array itself. 
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FIG. 6. Plot of the real part of dS2i/dH vs H for sample 
fla at an excitation frequency of 9.55 GHz, as measured by 
cavity FMR. The high-field mode seen in Figure[5]is preserved 
in the uniform excitation field of the FMR cavity. No units are 
included on the vertical axes because the scale is dependent 
on the particular tuning of the cavity. 



The relationship between the resonance field H^^s of 
both modes and the relative lattice-applied field angle (p 
was measured for sample f3b*. Neither mode showed any 
angular variation beyond the the 15 and 20 Oe scatter for 
the fundamental and high-field modes, respectively. This 
scatter was caused by a combination of small movements 
of the sample in the applied magnetic field when varying 
(j) and uncertainty in the baseline of the Lorentzian fit 
used to extract the values of iJres- This measurement 
restricts the angular variation of iJres to be below these 
values of uncertainty. Additionally, the amplitudes of 
both modes were essentially independent of cf). 

In previous studies by other authors on square arrays 
of sub-micron discs of similar diameter-to-pitch ratio, 
strong dependencies on the relative array-field angle 
of the resonance field have i?rcs have been measured for 
both fundamental and high-field modes.— ^^i^i Given the 
higher symmetry of the trigonal system, and the large 
diameter of the discs with respect to the pitch, it is not 
entirely surprising that such strong dependence was not 
observed in the fundamental mode of the arrays fabri- 
cated for this study. 

The / vs iJrcs data for film fla is shown in Figure [T] 
this data was representative of equivalent measurements 
of the other four samples. Resonances below 7 GHz were 
not recorded, so as to ensure that the resonance field 
was always high enough for the discs to be tangentially 
magnetically saturated. The data for each of the two 
modes for each sample was fit with Equation [1] under 
the assumption Nz = N^ , using Ny — Nz as a fit param- 
eter. The results of this process are tabulated in Table 
im There was no clear variation with (/)' in Ny — N^ for 
either mode across the samples. 

Significant differences between the resonance frequen- 
cies of close-packed and well-isolated discs have been re- 
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FIG. 7. Plot of resonance frequency / against resonance field 
Hres for sample fla (squares; circles - high-field mode) and 
flc (diamonds), a sample of the continuous film from which 
fla was patterned, as measured by VNA-FMR. No resonances 
were observed below the saturating field of the dots. 



Film 


Ny - TV^ fundamental Ny - iV^ high-field 


f3b* 


0.958 ± 0.001 


0.740 ± 0.002 


f3b 


0.934 ± 0.001 


0.716 ±0.001 


f3a 


0.939 ± 0.001 


0.720 ± 0.002 


f2a 


0.947 ± 0.001 


0.724 ± 0.001 


fla 


0.937 ± 0.001 


0.709 ± 0.002 



TABLE II. Table showing fit parameters Ny for the funda- 
mental and high-field modes, as extracted by fitting Equation 

m 



ported for square arrays in one other paper, -'^^ This effect 
was not observed in this study for the case where ar- 
ray ordering-rather than array pitch-was varied: there 
appears to be no correlation for this sample series be- 
tween array order cf)' and fitted demagnetizing parameter 
N„-N,. 



C. FMR linewidth 

The full- width half-maximum linewidth values AiJ ex- 
tracted from Lorentzian fits to the FMR modes are shown 
in Figure |51 Each plot shows the /S.H vs / data for the 
unpatterned parent film, and both the fundamental and 
high-field modes for the patterned films. The uncertain- 
ties in these linewidth values are not shown on the graphs 
for clarity: for the fundamental mode these uncertainties 
were on the order of ±2%, and for the high-field mode 
±20%. In both cases the primary source was uncertainty 
in the baseline of the Lorentzian fit; the uncertainty was 
more severe for the high-field mode because its smaller 
amplitude and proximity to the fundamental mode made 
the value of the linewidth highly dependent on the num- 



ber of points on either side of the peak used in the fit. It 
is clear from examination that the hnewidths of both the 
fundamental and high-field modes increase with increas- 
ing array disorder. 



Film 


0'[° mm-i] 


a X 10* 


A/Jo [Oe] 


Ai/o, a [Oc] 


f3c 


- 


9.28 ±0.16 


0.25 ±1.20 


- 


f3b* 


2.6 ±0.4 


8.61 ±0.16 


8.24 ±1.26 


2.31 ±3.10 


f3b 


6.0 ±0.8 


10.70 ±0.18 


1.02 ±1.38 


13.68 ± 3.33 


fSa 


9.4 ±1.1 


10.84 ±0.19 


3.54 ±1.48 


17.89 ±3.41 


f2c 


- 


7.79 ±0.14 


2.48 ± 1.06 


- 


f2a 


11.3 ±1.7 


11.00 ±0.19 


3.27 ±1.50 


32.21 ±3.25 


flc 


- 


7.77 ±0.14 


2.85 ± 1.06 


- 


fla 


19.9 ±2.1 


9.55 ±0.21 


19.90 ± 1.64 


36.08 ± 3.33 



TABLE III. Table showing fit data for each sample to Equa- 
tion [5] for the fundamental mode data, for the case where a is 
allowed to vary for patterned films and for the case when it is 
held to the value Ocnts of the corresponding continuous film. 



Film <^'[° mm-i] 


Q X 10* Affo [Oe] 


f3b* 2.6 ± 0.4 


5.00 ±1.14 14.74 ±9.23 


f3b 6.0 ± 0.8 


6.21 ±1.06 0.95 ±8.21 


f3a 9.4 ±1.1 


6.40 ±1.11 -0.92 ±8.36 


f2a 11.3 ±1.7 


6.10 ±1.12 2.80 ±8.70 


fla 19.9 ±2.1 


8.45 ±1.61 6.68 ±12.37 



TABLE IV. Table showing fit data for each sample to Equa- 
tion[2]for the high-field mode data. The large uncertainties of 
the fit parameters are the direct result of an estimated 20% 
error on high-field mode linewidth data. 

In the context of studies of spin-wave mode broaden- 
ing, Shaw et al. have pointed out the importance of sep- 
arating the effects of intrinsic damping from inhomoge- 
neous damping^i^ in the total linewidth A_ff, as defined 
by the equation: 



AiJ = Ai?o±4™/(7A^o/) (SI). 



(2) 



Here / is the frequency of precession and 7 the gy- 
romagnetic ratio, a is the intrinsic damping parame- 
ter in the Landau-Lifshitz-Gilbert equatiouj^i^ and can 
be thought of as a 'viscous' damping of energy to the 
latticed 

AiJo is a term representing inhomogeneous broaden- 
ing: an example of such a broadening might be the 
decay of uniform motion into spin waves having non- 
zero wavevectors, which might themselves decay to the 
lattice.?^- The IS.H vs / data represented in Figure |S](a)- 
(d) were fit with equation^ The extracted fit parameters 
a and AiJp are given in Table IIIII for the fundamental 
modes and Table ITVl for the high-field modes. 

Previous studies by other authors of FMR linewidth 
broadening in disc-geometry nanostructures?^^'— have 



produced some conflicting results as to whether pattern- 
ing affects the value of the intrinsic damping parameter, 
a. In the study presented here, a values did appear to 
change outside the bounds of experimental uncertainty. 
For the fundamental mode as measured in this study, a 
was larger for the patterned samples than for their cor- 
responding parent continuous films; the exception was 
the most ordered sample, for which a was slightly lower. 
Conversely, a for the high-field mode was lower than or 
within the range of experimental error of the correspond- 
ing value for the fundamental mode. Although the values 
of a for patterned films appeared to be different than for 
the corresponding parent films, there was no clear vari- 
ation of a for either mode with the parameter (j)' which 
distinguished the films. 

It is difficult to interpret these a values because of 
the lack of a complete theoretical description in the lit- 
erature of spin wave modes in closely packed and po- 
tentially dipole coupled tangentially magnetized disc ar- 
rays. This lack is a consequence of the difficulties caused 
by the absence of cylindrical symmetry in the plane of 
the discs,— >^ and complicated for this experimental sys- 
tem by the lack of perfect ordering. For these reasons, 
it is difficult to know that Equation [2] holds in this sys- 
tem. In particular, one can speculate that the inhomoge- 
nous broadening A7?o may have a frequency- or field- 
dependence which varies with array ordering, in anal- 
ogy to or even contributed to by the field dependence 
of extrinsic damping related to two-magnon scattering 
processes. ^^ Such dependency might go some way to ex- 
plaining the apparent differences in a between patterned 
and unpatterned films. If Ai/o is a function of field (or 
frequency) , then the steeper gradients in /S.H seen in Fig- 
ure [8] might not require the 'real' intrinsic damping pa- 
rameter to differ between patterned and parent films. 

Since the origin of these differing gradients was not 
clear, a second set of fits to Equation [2] was performed, 
with the slopes of the patterned film data set to be equal 
to those of the parent continuous films: a = ac- The 
difference A_ffo, a — Ai?o, c between this 'fixed intrinsic 
damping inhomogeneous broadening' value AiJo, a and 
the continuous film inhomogeneous broadening, ATJq, c 
is equivalent to the average of the difference between 
linewidths /S.H of the patterned and unpatterned films 
over the range of frequencies measured. It is a conve- 
nient experimental index of the inhomogeneous broad- 
ening that takes into account possible changes to a in 
the frequency domain of interest, rather than a quanti- 
tative measure of the 'true' inhomogeneous broadening 
as defined by Equation [2] A plot of this difference of in- 
homogeneous broadening, Aifo_ a ~ AiJo, c, vs ordering 
parameter 0' is shown in Figure |9l The increase of in- 
homogeneous broadening with increasing array disorder 
can be clearly seen. 

The difficulties encountered in fitting the Ai/o vs / 
data were even more pronounced for the high-field mode. 
Since this mode is unique to the patterned geometry, 
there was no analogue from the parent continuous films 
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FIG. 8. Plot of full-width-half maximum linewidth AH against resonance frequency / as measured by VNA-FMR for samples 
(a) f3b*; (b) f3b (open triangles: fundamental mode; open pentagons: high-field mode); (c) f2a; and (d) fla. The filled diamonds 
are the data points for the parent continuous films. For all samples except f3b, the open squares denote the fundamental mode, 
and the open circles the high-field mode of the patterned films listed; the key for fSb has already been listed in parentheses. 



with which to compare data. Additionally, there was 
a high degree of experiniental uncertainty in high-field 
mode AH values, which led to high uncertainties in a. 
This in turn produced AHq uncertainties that were com- 
parable to or greater than the AHq values themselves. 
Despite these uncertainties, the data in Figure |S] suggest 
that the linewidth for this mode also increases with in- 
creasing (f)' in the frequency range measured. 

Disc uniformity has been maintained as far as exper- 
imentally practicable, as shown by consistency of hys- 
teretic measurements and AFM images. Differences in 
linewidth and saturation magnetization of parent films 
have been accounted for in the comparison between pat- 
terned films. The variation in broadening with changing 
array ordering suggests that the discs are dipole coupled 
together, and that the degree of coupling affects the FMR 
linewidth of both modes. That sample f3b* shows a more 
narrow linewidth than the patterned film f3b from which 
it was cut supports this conjecture. It might have been 
speculated that differently arranged array domains in the 
more disordered samples, interacting strongly internally 
but weakly with other domains, were producing differ- 
ent resonance fields due to angular dependence of -ffrcs- 



However, the lack of angular dependence of -ffros in either 
mode in sample f3b* suggests that the broadening is not 
due to the sum of differing resonance fields of different 
array domains, but rather to some less local effect. 

In the only previous study of the FMR characteristics 
of trigonal arrays fabricated by a self-assembly process, 
it was not clear whether the anomalously large linewidth 
recorded was due to the fabrication process, a distribu- 
tion of disc parameters, or the quality of the parent con- 
tinuous film.^^ The results of the present study indicate 
that there is no dramatic impact on linewidth due to pat- 
terning in such a process: even for the most disordered 
array, fla, there was less than a factor-of-two broaden- 
ing; for the case of near-perfect order represented by 
sample f3b*, there was essentially no broadening. The 
broad linewidth seen in Reference [ll| is therefore likely 
to have been caused by a broad linewidth in the parent 
film rather than the patterning process itself. The nar- 
rowness of the 'waist' of the vortex hysteresis loops and 
high magnetization of the patterned films studied here 
are consistent with this conclusion. The utility of the 
water-surface self-assembly lithographic process for pat- 
terning disc arrays has thus been partially vindicated: 
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FIG. 9. Plot of AHo, a - A.H'o.c vs (j)' for the fundamental 
mode data of patterned films f3b*, f3b, f3a, f2a, and fla. 
The data plotted here corresponds to the appropriate data 
presented in tabular format in Table UTIl 



despite the reasonable criticisms that the self-assembly 
process does not preserve long range order, the impact of 
variations in that order on the FMR linewidth have been 
shown to be relatively minor. 



IV. CONCLUSIONS 

A series of arrays of sub-micron Permalloy discs was 
fabricated by a water-surface self-assembly lithographic 
process. The geometry, static magnetic behavior, and 
ordering of the disc arrays were measured. The arrays 
were characterized by vector network analyzer ferromag- 
netic resonance. Two modes were present in the ferro- 
magnetic resonance spectra: a fundamental mode and a 
higher field, lower frequency mode. The variation of the 
resonance field with relative lattice-to-static field angle 
was shown to be negligible for both modes. The rela- 
tionship between resonance frequency and resonance field 
was found to be uncorrelated to increasing array disorder. 
The linewidth vs frequency data for these arrays revealed 
that the self-assembly fabrication process did not dra- 
matically increase the ferromagnetic resonance linewidth, 
clarifying questions from a previous study. The intrinsic 
damping parameter for the fundamental mode appeared 
to be higher for the patterned films than in the parent 
continuous films, but it was not clear whether this was 
a result of field- or frequency-dependencies of the inho- 
mogeneous damping not accounted for in the fit. The 
linewidth of both the fundamental and high-field modes 
was found to increase with increasing array disorder. 
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